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In  this  work,  we  combine  electrochemical  measurements,  such  as  the  rotating  disk  electrode  (RDE)  and  the 
rotating  ring-disk  electrode  (RRDE)  techniques,  and  density  functional  theory  (DFT)  calculations  to  elucidate 
the  mechanisms  of  the  oxygen  reduction  reaction  (ORR)  on  carbon-supported  Fe-phthalocyanine  (FePc/C) 
and  Co-phthalocyanine  (CoPc/C)  catalysts  in  0.1  M  NaOH  solutions.  The  onset  potential  for  ORR  on  FePc/C 
catalyst  is  found  to  be  around  0.05  V  vs.  Hg/HgO  in  0.1  M  NaOH  solutions,  which  is  100  mV  more  positive 
than  that  on  CoPc/C.  RDE  and  RRDE  measurements  show  that  the  ORR  mechanism  is  via  a  4e“  pathway  on 
the  FePc/C  while  it  is  through  a  2e“  pathway  on  the  CoPc/C  catalyst.  The  catalyst  stability  tests  reveal  that 
FePc  is  much  less  stable  than  CoPc  under  fuel  cell  cathode  working  conditions.  Moreover,  DFT  calculations 
were  performed  to  study  the  adsorption  of  O2,  H20,  OH,  HOOH,  and  H200  molecules  on  FePc  and  CoPc 
molecule  catalysts.  We  conclude  the  following  from  our  theoretical  and  experimental  results  for  the  ORR  on 
FePc/C  and  CoPc/C  catalysts:  (1)  the  lower  the  02  adsorption  energy,  the  higher  the  kinetics  of  the  ORR  will 
be;  (2)  the  ORR  pathways,  2e“  or  4e“,  are  mainly  determined  by  the  H202  adsorption  configurations;  and  (3) 
OH  adsorption  on  the  catalysts  is  considered  to  be  an  important  factor  to  affect  the  catalyst  stability. 


I.  Introduction 

Recently  there  has  been  a  growing  worldwide  interest  in 
developing  solid  alkaline  fuel  cells  (SAFCs),  which  have  various 
applications  as  low-cost  yet  high-performance  power  sources.1-10 
Compared  to  liquid  alkaline  fuel  cells  (AFCs)  and  proton 
exchange  membrane  fuel  cells  (PEMFCs),  SAFCs  possess 
several  important  merits,  including  (1)  facile  kinetics  at  the 
cathode  and  the  anode,  (2)  minimal  methanol/ethanol  crossover 
effect,  (3)  simplified  water  management,  (4)  no  precipitated 
carbonate,  (5)  no  electrolyte  weeping,  and  (6)  a  wide  choice  of 
materials  for  current  supports,  bipolar  plates,  and  fuel  cells.  To 
advance  the  technology  of  SAFCs,  we  have  been  developing 
novel  alkaline  anion  exchange  membranes1112  and  investigating 
various  transition  metals13,14  as  well  as  transition  metals  of  N4- 
macrocycles  (M-N4-macrocycles)  as  electrocatalysts  to  replace 
Pt  for  the  oxygen  reduction  reaction  (ORR)  in  alkaline  media. 

M-N4-macrocycles,  such  as  porphyrins  and  phthalocyanines, 
have  been  considered  as  promising  catalysts  for  the  ORR15-28 
in  alkaline  and  acid  media  for  several  decades.  It  has  been 
established  that  (1)  the  type  of  central  metal  atom  determina- 
tively  influences  their  catalytic  activity,  (2)  their  polymeric  or 
heat-treated  M-N4  forms  could  lead  to  better  chemical  and 
thermal  stabilities,  higher  electronic  conductivities,  and  higher 
02  activities23-25  than  their  monomeric  species,  and  (3)  the  pH 
values  of  electrolytes  affect  the  activities  and  stabilities  of  the 
catalysts.16  Despite  much  effort  in  developing  M-N4-macrocycles 
as  fuel  cell  electrocatalysts,  so  far,  no  macrocycle  molecule 


*  To  whom  correspondence  should  be  addressed.  E-mail:  rochen@iupui.edu 
(R.C.)  or  wang83@iupui.edu  (G.W.). 

1  Richard  G.  Lugar  Center  for  Renewable  Energy,  Indiana  University 
Purdue  University  Indianapolis. 

*  U.S.  Army  Research  Laboratory. 

§  Department  of  Mechanical  Engineering,  Indiana  University  Purdue 
University  Indianapolis. 


catalysts  have  demonstrated  acceptable  ORR  catalytic  activity 
and  stability  as  compared  to  Pt-based  catalysts.  At  the  moment, 
lacking  a  fundamental  understanding  of  the  catalytic  activity 
and  stability  of  these  catalyst  molecules  hinders  further  develop¬ 
ment  for  their  practical  applications.  In  this  paper,  we  combine 
electrochemical  measurements,  such  as  the  rotating  disk  elec¬ 
trode  (RDE)  and  the  rotating  ring-disk  electrode  (RRDE) 
techniques,  and  density  functional  theory  (DFT)  calculations 
to  elucidate  the  mechanisms  of  the  ORR  on  carbon-supported 
Fe-phthalocyanine  (FePc/C)  and  Co-phthalocyanine  (CoPc/C) 
in  alkaline  media.  The  synergic  combination  of  theory  and 
experiments  enables  us  to  better  interpret  the  ORR  overpotential, 
the  reaction  pathways,  and  the  stability  of  these  catalysts. 

The  electrochemical  reduction  of  02  is  a  multielectron 
reaction  that  has  two  main  possible  pathways:  one  involving 
the  transfer  of  two  electrons  to  produce  H202,  and  the  other,  a 
direct  four-electron  pathway  to  produce  water.  In  alkaline  media, 
the  2e-  pathway  can  be  written25  as: 

02  +  H20  +  2e"  —  H02  +  OH 

E0  =  -0.230  V  vs.  Hg/HgO  (0. 1  M  NaOH)  (1) 

The  direct  4e~  pathway  can  be  written  as: 

02  +  2H20  +  4e~  —  40H 

E0  =  0.236  V  vs.  Hg/HgO  (0. 1  M  NaOH)  (2) 

To  obtain  maximum  energy  capacity,  it  is  highly  desirable 
to  reduce  02  via  the  4c  pathway.  Precious  Pt  and  its  alloys 
are  known  to  promote  the  4e  ORR  pathway.  For  Pt  and  other 
metallic  catalysts,  DFT  calculations  have  been  successfully 
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applied  to  elucidate  the  origin  of  the  overpotential  of  the  ORR,29 
the  mechanism  of  the  ORR  in  acids,30  and  the  alloying  effect 
on  the  catalytic  activity  for  the  ORR.31  However,  a  similar  study 
is  scarce  for  other  catalyst  systems  like  M-N4-macrocycles. 
Recently  we  reported  our  work32  studying  the  adsorption  of 
oxygen  molecules  on  FePc  and  CoPc  molecules  using  the  DFT 
method.  We  found  that  there  are  multiple  stable  adsorption 
configurations  for  O2  adsorbed  on  FePc  and  CoPc  molecules. 
Extending  our  prior  work,  here  we  performed  detailed  DFT 
calculations  for  O2,  HzO,  OH,  and  H2O2  molecule  adsorption 
on  the  surfaces  of  FePc  and  CoPc  in  order  to  gain  insights  into 
the  ORR  mechanism  on  those  catalyst  molecules. 

II.  Experimental  and  Theoretical  Methods 

11.1.  Materials.  The  FePc  and  CoPc  were  purchased  from 
Alfa  Aesar  and  used  as  received.  To  disperse  FePc  or  CoPc  on 
a  high  surface  area  carbon  substrate  (Vulcan  XC-72,  CABOT 
Corp.,  Boston,  MA,  USA),  86.340  mg  of  FePc  or  CoPc  was 
first  weighted,  dissolved  in  concentrated  H2SO4,  and  sonicated 
for  30  min.  After  complete  dissolution  of  the  FePc  or  CoPc  in 
the  H2SO4  solution,  0.2  g  of  Vulcan  XC-72  carbon  supports 
was  added  into  the  solution  which  was  ultrasonically  homo¬ 
genized  to  form  a  suspension  of  FePc  or  CoPc  dispersed  on 
carbon  supports.  Then  the  suspension  was  poured  into  cold 
deionized  (DI)  water,  filtered,  washed  with  DI  water  again  to 
remove  H2S04  residuals,  and  finally  dried  at  70  °C  in  an  oven 
overnight  under  ambient  air  conditions.  It  is  worth  mentioning 
that  organic  DMF  (dimethylformamide)  solvent  can  replace 
concentrated  H2S04  in  the  above  procedures  to  dissolve  FePc 
or  CoPc  and  lead  to  similar  electrochemical  measurement  results 
for  the  prepared  catalysts.  In  this  work,  we  found  that  the  carbon- 
supported  FePc  or  CoPC  catalysts  prepared  with  the  acid  solvent 
have  slightly  better  performance  than  those  prepared  with  DMF 
solvent.  Consequently,  we  here  only  present  the  results  for  the 
samples  prepared  with  H2SO4. 

The  carbon-supported  electrocatalyst  Pt/Vulcan  XC-72  (20 
wt  %  metal  on  carbon)  was  purchased  from  BASF  and  used  as 
received.  Sodium  hydroxide  (J.T.  Baker)  was  purchased  with  a 
certified  ACS  reagent  grade  and  used  as  received.  O2  gas  was 
purchased  from  Praxair.  All  aqueous  suspensions  used  in  this 
study  were  prepared  with  Millipore  ultrapure  water  (18.2 
M£2*cm). 

FePc/C  or  CoPc/C  ink  solutions  were  prepared  by  mixing  3 
mg  of  the  FePc/C  or  CoPc/C  catalyst  with  1  mL  of  ethanol,  1 
mL  of  water,  and  33  /tL  of  0.5%  Nafion  in  an  ultrasonic  bath 
for  1  h.  Comparative  Pt/C  (20%  weight)  ink  was  prepared  by 
mixing  5  mg  of  Pt/C  catalyst  with  1.25  mL  of  ethanol,  1.25 
mL  of  water,  and  52.6  /tL  of  0.5%  Nafion  in  an  ultrasonic  bath 
for  1  h. 

11.2.  Preparation  of  Electrodes.  A  grassy  carbon  (GC)  disk 
electrode  (Pine  Instruments)  with  a  5.0  mm  diameter  was  used 
as  the  working  electrode.  Before  any  measurements,  the 
electrode  surface  was  polished  with  an  aqueous  suspension  of 
ALO3  (0.05  Mtn)  and  rinsed  with  DI  water.  The  GC  electrode 
was  then  placed  in  an  ultrasonic  cleaning  bath  for  5  min  and 
rinsed  again  with  DI  water.  To  absorb  the  molecule  catalysts 
on  it,  the  GC  disk  electrode  was  immersed  into  a  dimethylfor¬ 
mamide  (DMF)  solution  containing  1.0  x  10  3  M  FePc  or  CoPc 
and  rotated  at  500  rpm  for  25  min.  Then  the  electrode  was  rinsed 
with  DI  water  and  introduced  into  the  electrochemical  cell  for 
the  ORR  stability  test. 

The  working  electrodes  of  FePc/C,  CoPc/C,  or  Pt/C  were 
made  by  dropping  40  fiL  of  the  corresponding  ink  solution  on 
the  surface  of  either  the  GC  rotating  disk  electrode  or  the  rotating 


ring-disk  electrode.  Then  the  deposited  drop  of  ink  was  dried 
in  an  air-flow  environment  at  room  temperature  to  evaporate 
the  solvent.  Thus,  an  ultrathin  film  of  FePc/C,  CoPc/C,  or  Pt/C 
catalyst  was  formed  on  the  GC  surface.  In  this  work,  the  three 
catalyst  films  have  a  nearly  equal  thickness,  which  is  about  1.70 
um.  On  the  basis  of  their  densities,  the  metal  loading  on  the 
disk  electrode  was  estimated  to  be  about  8.7  Mg/cm2  for  Fe,  9.1 
fig/cm2  for  Co,  and  83  /tg/cm2  for  Pt.  It  should  be  noted  that 
the  metal  loadings  are  estimated  here  with  the  assumption  that 
the  catalyst  molecules  (FePc  or  CoPc)  are  completely  dissolved 
in  the  concentrated  H2SO4  acid  solution  and  distributed  uni¬ 
formly  in  the  ink  when  preparing  the  thin  films.  In  this  work, 
we  did  not  conduct  further  measurements  on  the  exact  composi¬ 
tion  of  the  catalysts  during  the  electrochemical  tests. 

11.3.  Electrochemical  Measurements.  The  electrochemical 
measurement  setup  consists  of  either  a  computer-controlled 
Solartron  1287  electrochemical  interface  or  Pine  bipotentiostat, 
a  radiometer  speed  control  unit  from  Pine  Instrument  Company 
(MSRX  Speed  Control),  and  a  rotating  ring-disk  electrode 
radiometer  (RRDE,  glassy  carbon  with  a  diameter  of  5.0  mm 
as  the  disk  and  Pt  as  the  ring).  The  RDE  and  the  RRDE 
experiments  were  conducted  in  an  02-purged  0.1  M  NaOH 
solution  using  a  conventional  standard  three-electrode  electro¬ 
chemical  cell  with  a  jacket.  A  Pt  wire  counter  electrode  with  a 
large  surface  area  and  a  Hg/HgO  reference  electrode  connected 
via  a  salt  bridge  were  used.  All  the  potentials  in  this  work  refer 
to  Hg/HgO  in  a  0.1  M  NaOH  solution,  denoted  as  Hg/HgO 
(about  0.164  V  vs.  NHE). 

H2O2  production  in  02-saturated  0.1  M  NaOH  electrolytes 
was  monitored  in  a  RRDE  configuration,  using  a  polycrystalline 
Pt  ring  biased  at  0.3  V  vs.  Hg/HgO.  The  ring  current  (7r)  was 
recorded  simultaneously  with  the  disk  current  (/d).  Collection 
efficiency  (V)  of  the  ring  electrode  was  calibrated  by  K3Fe(CN)6 
redox  reaction  in  an  Ar-saturated  0.1  M  NaOH  solution.  The 
value  of  the  collection  efficiency  (V  =  7r/7d)  determined  is  0.23 
for  both  the  FePc/C  and  the  CoPc/C,  and  0.26  for  the  Pt/C 
electrodes.  The  fractional  yields  of  H2O2  in  the  ORR  were 
calculated  from  the  RRDE  experiments  as  XH2o2  =  (27r//V)/(/d 
+  2/r//V).  The  mass  transport  normalized  kinetic  ORR  currents 
at  the  disk  electrode  were  calculated  as  7k  =  7lim//(7Um  —  7),  with 
limiting  current  (7iim)  taken  as  the  ORR  current  at  potential  —0.4 
or  —0.5  or  —0.6  V  vs.  Hg/HgO,  respectively. 

11.4.  Computational  Methods.  DFT  calculations  were  per¬ 
formed  with  the  Materials  Studio  DMol3  module  from  Accel- 
rys.33  The  generalized  gradient  approximation  of  Perdew  and 
Wang34  was  used  for  exchange  and  correlation  and  the  double 
numeric  basis  with  polarization  functions  was  used  as  the  atomic 
basis  set.  Density  functional  semicore  pseudopotentials,  which 
were  generated  by  fitting  all  electron  relativistic  DFT  results, 
were  used  to  reduce  the  computation  cost.  Spin-unrestricted 
wave  functions  and  Fermi  orbital  occupations  were  used  in  our 
calculations.  The  energy  convergence  for  geometry  optimization 
was  set  to  1  x  1CT5  eV. 

III.  Results 

III.l.  RDE  and  RRDE  Measurements.  Several  important 
reaction  parameters,  such  as  the  number  of  electrons  involved 
in  the  ORR,  Tafel  slopes,  and  the  kinetic  rate  constant  for  the 
ORR,  can  be  obtained  from  RDE  measurements.  Figure  1  shows 
the  rotating  disk  polarization  curves  at  various  rotation  rates 
for  oxygen  reduction  on  a  GC  RDE  electrode  coated  with  the 
FePc/C  (Figure  la),  the  CoPc/C  (Figure  lb),  or  the  20  wt  % 
Pt/C  (Figure  lc)  catalyst  in  a  0.1  M  NaOH  solution  saturated 
with  oxygen.  The  onset  ORR  potential  on  the  FePc/C  electrode 
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electrode  in  a  1  M  KOH  solution,  and  Tanaka  et  al.37  showed 
the  onset  potential  of  02  reduction  to  be  0.050  V  vs.  Hg/HgO 
in  a  1  M  NaOH  solution.  Jiang  et  al.38  reported  that  the  onset 
potential  for  02  reduction  on  the  Pt/C  is  around  0.050  V  vs. 
Hg/HgO  in  a  0.1  M  NaOH  solution,  corresponding  to  an 
overpotential  of  0.190  V.  M.  Issacs  et  al.39  showed  the  onset 
ORR  potential  on  a  CoPc  electrode  around  —0.150  V  vs.  SCE 
in  a  0.1  M  NaOH  solution. 

Unlike  in  0. 1  M  H2S04  solution,27  the  polarization  curves  of 
the  FePc/C  electrode  in  0.1  M  NaOH  solution  (see.  Figure  la) 
display  current  plateaus  in  the  high-polarization  range.  The  mass 
transport  limiting  ORR  current  was  reached  at  the  potential  less 
than  —0.1  V  vs.  Hg/HgO  on  the  FePc/C  and  less  than  —0.2  V 
vs.  Hg/HgO  observed  on  the  CoPc/C.  The  mixed  mass  transport 
and  kinetic  controlled  region  was  reached  at  the  potentials  from 
—0.1  to  0.05  V  on  the  FePc/C,  from  —0.2  to  —0.05  V  on  the 
CoPc/C,  and  from  —0.2  to  0.05  V  on  the  Pt/C.  The  half-wave 
potentials  were  around  —0.050  V  on  the  FePc/C,  —0.100  V  on 
the  Pt/C,  and  —0.120  V  on  the  CoPc/C.  Our  results  here  also 
agree  well  with  what  had  been  reported.35,38,39  To  examine  the 
influence  of  Vulcan  carbon  on  the  catalytic  performance  of  the 
catalysts,  RDE  experiments  were  also  conducted  on  the  elec¬ 
trodes  with  noncarbon-supported  CoPc  and  FePc  catalysts  as 
well  as  with  just  Vulcan  carbon.  Comparing  the  curves  for  the 
same  catalysts  with  and  without  carbon  supports,  we  found  that 
Vulcan-carbon  contributes  little  toward  the  ORR  current  for  the 
FePc/C  catalysts  but  contributes  appreciably  to  the  ORR  current 
for  the  CoPc/C  catalysts. 

The  Koutecky— Fevich  equation40  states  that 

7lim  =  0.62  nFD2,3v~V6c0a)112  (3) 
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Figure  1.  RDE  polarization  curves  at  different  rotation  rates  for  oxygen 
reduction  on  GC  with  predeposited  (a)  FePc/C,  (b)  CoPc/C,  and  (c) 
Pt/C  in  a  0. 1  M  NaOH  solution  saturated  with  oxygen.  GC  electrode 
area  =  0.196  cm2,  scan  rate  =  10  mV  s-1,  and  room  temperature.  In 
the  figure,  the  potentials  are  given  with  reference  to  both  Hg/HgO 
electrode  and  normal  hydrogen  electrode  (NHE).  Note  that  the  computed 
potential  of  the  Hg/HgO  electrode  versus  an  NHE  at  25  °C  and  1  atm 
is  0.164  V  for  O.f  M  NaOH  solution. 

starts  at  the  same  potential  (0.050  V  vs.  Hg/HgO)  as  that  on 
the  Pt/C  electrode  and  shifts  more  than  100  mV  positively  as 
compared  to  that  on  the  CoPc/C  electrode.  These  results  are  in 
good  agreement  with  previous  ones  reported  in  the  literature.35-39 
For  example,  Behret  et  al.35  reported  the  onset  ORR  potential 
to  be  0.9  V  vs.  dynamic  hydrogen  electrode  (DHE,  —960  mV 
vs.  Hg/HgO  in  5  M  KOH36)  measured  with  a  rotating  FePc  disk 


where  /um  is  the  limiting  current  density,  n  is  the  number  of 
electrons  transferred  per  oxygen  molecule,  F  is  the  Faraday 
constant  (96500  C  mol-1),  D  is  the  02  diffusion  coefficient  (1.9 
x  10-5  cm2  s-1)  in  0.1  M  NaOH,  v  is  the  kinematic  viscosity 
(1.1  x  10-2  cm2  s-1),  Co  is  the  concentration  of  oxygen  (1.2  x 
10-6  mol  cm-3),  and  a>  is  the  rotation  rate  in  the  radian.16,39,41,42 
Figure  2  shows  the  Koutecky— Levich  plots  of  l//iim  vs.  llcum 
at  fixed  potentials  on  the  FePc/C  (Figure  2a),  the  CoPc/C  (Figure 
2b),  and  the  Pt/C  (Figure  2c)  electrodes  derived  from  the  data 
in  Figure  1.  From  the  slopes  of  the  Koutecky— Levich  plots, 
we  calculated  the  electron  exchange  number  ( n )  for  ORR  at 
the  fixed  potentials  of  —0.4,  —0.5,  and  —0.6  V  vs.  Hg/HgO 
(mass  transport  regime).  For  both  FePc/C  and  Pt/C,  n  is 
determined  to  be  3.8,  agreeing  well  with  what  has  been 
reported.16,38  This  result  indicates  that  the  02  reduction  catalyzed 
on  the  FePc/C  or  on  the  Pt/C  electrode  is  a  4e-  reduction  process 
leading  to  OH  as  the  main  product.  For  the  CoPc/C,  n  is  in  a 
range  of  1.4  to  2.4,  which  is  consistent  with  the  data  from 
another  study.39  This  result  suggests  that  the  ORR  catalyzed  on 
the  CoPc/C  electrode  is  a  2e-  reduction  process  leading  to  H202. 

To  further  verify  the  ORR  pathways  on  the  FePc/C  and  the 
CoPc/C  electrodes,  we  carried  out  RRDE  measurements  to 
monitor  the  formation  of  H202  during  the  ORR  process.  Figure 
3  gives  the  RRDE  polarization  curves  for  02  reduction  on  the 
FePc/C,  CoPc/C,  and  Pt/C  with  a  rotation  rate  of  2500  rpm  in 
0.1  M  NaOH  solution  saturated  with  oxygen.  The  inset  of  Figure 
3  gives  the  H202  yields  on  the  three  catalysts.  On  both  FePc/C 
and  Pt/C  electrodes,  no  significant  solution  phase  H202  was 
detected  and  thus  the  H202  yield  was  negligible.  For  the  CoPc/C 
electrode,  a  significant  ring  current  was  detected  starting  at  the 
ORR  onset  potential  of  the  disk  electrode  and  up  to  50%  of  the 
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Figure  2.  Plot  of  I1  vs.  G)-1/2  for  oxygen  reduction  on  GC  with 
precoated  (a)  FePc/C,  (b)  CoPc/C,  and  (c)  Pt/C  catalysts  in  a  0.1  M 
NaOH  solution  saturated  with  oxygen  at  potentials  —0.4,  —0.5,  and 
-0.6  V  vs.  Hg/HgO. 


H2O2  yield  was  measured.  This  result  indicates  that  H2O2  is  a 
main  product  for  the  ORR  catalyzed  by  the  CoPc/C  but  not  for 
ORR  catalyzed  by  the  FePc/C.  Referencing  eqs  1  and  2,  the 
RRDE  results  confirm  our  results  calculated  from  the  Levich 
equation  that  the  ORR  electron  exchange  number  is  3.8  (close 
to  4e  )  for  the  ORR  on  the  FePc/C  and  Pt/C,  but  about  2  (2c  ) 
for  the  ORR  on  the  CoPc/C. 

Figure  4  presents  the  mass-corrected  Tafel  plots  of  log  4  (mA 
cm-2)  vs.  the  electrode  potential  E  for  the  ORR  on  the  FePc/C, 
the  CoPc/C,  and  the  Pt/C  electrodes  in  an  02-saturated  0.1  M 
NaOH  solution.  These  Tafel  curves  were  derived  from  the 
polarization  curves  of  Figure  1  with  a  rotation  rate  of  2500  rpm. 
The  Tafel  slopes  (listed  in  Table  1)  for  the  FePc/C  and  Pt/C 
catalysts  are  quite  similar  and  can  be  divided  into  two  parts. 


i.e.,  about  50  mV  dec-1  at  the  lower  overpotential  region  (where 
E  >  0  V  vs.  Hg/HgO)  and  about  120  mV  dec-1  at  the  higher 
overpotential  region  (where  E  <  — 100  mV  vs.  Hg/HgO).  These 
results  correspond  well  with  those  reported  in  the  literature16,38 
and  imply  that  the  rate-determining  step  of  the  ORR  is  similar 
for  both  the  FePc/C  and  Pt/C  catalysts.  The  Tafel  slope  of  the 
CoPc/C  catalyst  is  around  120  mV  dec-1  at  the  higher 
overpotential  region  (where  E  <  —150  mV  vs.  Hg/HgO), 
corresponding  well  with  what  were  reported  by  M.  Issacs  et 
al.39  At  a  typical  alkaline  fuel  cell  cathode  working  potential 
(—0.05  V  vs  Hg/HgO,  equivalent  to  an  overpotential  of  0.280 
V),  the  ORR  kinetic  current  on  the  FePc/C  is  10  mA  cm-2, 
which  is  about  3  times  higher  than  that  of  the  Pt/C  electrode 
(3.1  mA  cm-2).  In  a  pure  oxygen  and  0.1  M  NaOH  electrolyte 
environment,  the  FePc/C  displays  better  performance  than  both 
the  20  wt  %  Pt/C  and  the  CoPc/C. 

A  rapid  decline  of  the  catalytic  activity  for  ORR  on  both 
unsupported  and  carbon-supported  M-N4-macrocycles  electro¬ 
catalysts  has  been  reported.25,27, 43-45  To  examine  the  stability 
of  FePc  and  CoPc  as  electrocatalysts  for  ORR,  we  performed 
potentiostatic  measurements  with  a  potential  set  at  —0.100  V 
vs.  Hg/HgO  on  the  FePc/C  and  the  CoPc/C  electrodes  with  a 
rotating  rate  of  900  rpm  in  02-saturated  0.1  M  NaOH  solutions. 
The  ORR  current  density  at  the  held  potential  was  recorded  for 
6000  s  and  shown  in  Figure  5a.  The  ORR  current  density  on 
the  FePc/C  electrode  was  found  to  be  more  than  an  order  of 
magnitude  higher  than  that  on  the  CoPc/C,  but  declined  quickly 
with  test  time  at  a  rate  of  0.1  uA  cm-2  s  '.  In  contrast,  the 
ORR  current  density  on  the  CoPc/C  electrode  remains  constant. 
It  infers  that  the  FePc/C  does  not  have  as  good  stability  as  the 
CoPc/C  catalyst  for  ORR  in  alkaline  solutions.  However,  the 
comparison  in  Figure  5a  between  the  FePc/C  and  CoPc/C  in 
their  stability  is  not  convincingly  conclusive  because  the  carbon- 
supported  FePc  and  CoPc  may  degrade  differently  under 
different  operation  currents  and  cause  the  observed  discrepancy. 
To  address  that  problem,  we  further  performed  RDE  cycling 
tests  on  the  FePc  and  CoPc  catalysts  without  carbon  supports 
adsorbed  to  a  GC  electrode  with  a  scanning  rate  at  10  mV/s 
and  a  rotating  rate  of  2500  rpm  in  oxygen-saturated  0.1  M  NaOH 
solutions.  As  shown  in  Figure  5b,  the  FePc  catalyst  without 
carbon  supports  also  displayed  a  significant  reduction  in  its  ORR 
current  density  with  each  continuing  cycle,  while  the  CoPc 
catalyst  showed  less  obvious  ORR  current  density  reductions. 
Hence,  our  results  show  clearly  that  the  FePc  with  or  without 
carbon  supports  has  less  ORR  electrocatalytic  stability  than  the 
CoPc  under  fuel  cell  operation  conditions. 

A  performance  recovery  experiment  was  conducted  on  the 
FePc  catalyst  preadsorbed  on  a  GC  RDE  surface.  Figure  6 
displays  the  results  from  this  experiment.  We  first  run  a  cyclic 
voltammogram  test  on  the  RDE  electrode  preadsorbed  FePc 
molecules  for  30  cycles  by  sweeping  potentials  from  —0.700 
to  0.100  V  vs.  Hg/HgO  in  an  oxygen-saturated  0.1  M  NaOH 
solution  with  scanning  rate  at  10  mV  s_1.  As  shown  in  Figure 
6,  the  catalytic  activity  of  the  FePc  catalyst  dropped  dramatically 
at  the  end  of  the  30th  cycle.  Then,  we  left  the  electrode  in  the 
same  alkaline  solution  for  12  h.  Further,  a  new  electrochemical 
measurement  on  the  same  electrode  was  performed  in  the  0.1 
M  NaOH  solution  saturated  with  oxygen.  A  noticeable  recovery 
in  the  catalytic  activity  of  the  FePc  catalyst  for  ORR  could  be 
clearly  observed  in  Figure  6.  However,  the  catalytic  performance 
of  the  FePc  could  not  be  fully  recovered. 

III.2.  DFT  Calculations.  According  to  eqs  1  and  2,  we  know 
that  02  and  H20  molecules  are  reactants,  OH  molecule  is  a 
product  via  both  the  4c  and  the  2e  paths,  while  HOOH 
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Figure  3.  RRDE  measurements  of  the  ORR  on  the  FePc/C,  CoPc/C,  and  Pt/C  electrodes  with  a  rotation  rate  of  2500  rpm  in  a  0.1  M  NaOH 
solution  during  positive  scanning.  Scan  rate  =10  mV  s-1.  In  the  figure,  the  potentials  are  given  with  reference  to  both  Hg/HgO  electrode  and 
normal  hydrogen  electrode  (NHE).  Note  that  the  computed  potential  of  the  Hg/HgO  electrode  versus  an  NHE  at  25  °C  and  1  atm  is  0.164  V  for 
0.1  M  NaOH  solution. 


log  [lk  (mA/cm2)] 


Figure  4.  Tafel  plots  of  log  4  vs.  E  (V)  for  the  ORR  on  the  FePc/C, 
CoPc/C,  and  Pt/C  electrodes  in  an  oxygen-saturated  0.1  M  NaOH 
solution  at  a  rotation  rate  of  2500  rpm.  Limiting  current  density  obtained 
at  potential  —0.5  V  vs.  Hg/HgO  was  used  to  calculate  4-  GC  electrode 
area  =  0.196  cm2,  room  temperature.  In  the  figure,  the  potentials  are 
given  with  reference  to  both  Hg/HgO  electrode  and  normal  hydrogen 
electrode  (NHE).  Note  that  the  computed  potential  of  the  Hg/HgO 
electrode  versus  an  NHE  at  25  °C  and  1  atm  is  0.164  V  for  0.1  M 
NaOH  solution. 

(hydrogen  peroxide)  and  H200  (oxywater)  are  the  two  isomers 
of  H202,  which  is  another  product  molecule  via  the  two-electron 
route.  To  elucidate  the  reaction  pathway  of  the  ORR  on  FePc 
and  CoPc  catalyst  molecules,  we  performed  DFT  calculations 
to  determine  the  structure  and  energy  of  02,  FFO,  OH,  HOOH, 
and  HsOO  molecule  adsorption  on  the  surfaces  of  the  FePc  and 
the  CoPc  molecules. 


TABLE  1:  Tafel  Slopes  (mV  dec1)  for  FePc/C,  CoPc/C, 
and  Pt/C  at  Different  Voltages 

catalysts  low  rj  (at  E  >  0  V)  high  (E  <  —0.100  V) 

FePc/C  42  120 

CoPc/C  127 

Pt/C  47  128 

The  optimized  structure  for  the  above  adsorption  systems  was 
obtained  through  structural  optimization  calculations  using  the 
DFT  method  implemented  in  Dmol3.  After  attaining  the 
optimized  structures,  we  calculated  the  adsorption  energy  for 
those  molecules  on  the  phthalocyanine  molecules.  The  adsorp¬ 
tion  energy  is  defined  as  the  energy  difference  between  the 
adsorption  and  the  isolated  systems.  Here,  the  energy  of  the 
isolated  system  refers  to  the  sum  of  the  energies  of  the  relaxed 
phthalocyanine  and  the  individual  adsorbate  molecules.  There¬ 
fore,  negative  adsorption  energy  indicates  that  the  adsorbate 
molecule  would  be  energetically  favorable  to  be  adducted  to 
the  surface  of  the  FePc  and  the  CoPc  catalysts,  while  positive 
adsorption  energy  means  that  the  adsorbate  molecule  would  not 
be  energetically  favorable  to  make  contact  with  the  FePc  and 
the  CoPc  catalysts.  Consequently,  the  negative  adsorption  energy 
of  the  02  molecule  on  the  surface  of  the  catalysts  is  a 
requirement  for  the  catalysts  to  promote  ORR. 

In  Table  2,  we  give  the  predicted  adsorption  energy  for  the 
most  stable  (namely,  the  lowest  energy)  adsorption  configuration 
of  those  molecules  on  FePc  and  CoPc.  Moreover,  we  plotted 
in  Figure  7  the  optimized  lowest  energy  molecular  structures 
of  02,  H20,  and  OH  adsorbed  on  FePc  or  CoPc.  It  can  be  noted 
that  02,  H20,  and  OH  molecules  assume  the  same  configuration 
as  their  most  stable  adsorption  on  both  FePc  and  CoPc 
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Figure  5.  The  stability  test  for  FePc  and  CoPc  precoated  on  a  GC 
electrode  at  a  rotating  rate  of  900  rpm  in  oxygen-saturated  0.1  M 
NaOH  at  room  temperature,  (a)  Potentiostatic  measurements  at  —0.1 

V  vs.  Hg/HgO  on  FePc/C  and  CoPc/C  at  900  rpm;  (b)  cycling  test 
on  FePc  and  CoPc  preadsorbed  on  a  GC  electrode  with  a  scanning 
rate  10  mV  s*1  at  a  rotating  rate  of  2500  rpm.  In  the  figure,  the 
potentials  are  given  with  reference  to  both  Hg/HgO  electrode  and 
normal  hydrogen  electrode  (NHE).  Note  that  the  computed  potential 
of  the  Hg/HgO  electrode  versus  an  NHE  at  25  °C  and  1  atm  is  0.164 

V  for  0.1  M  NaOH  solution. 

molecules.  In  the  most  stable  adsorption  of  02,  one  O  atom 
lies  right  above  the  central  metal  atom  over  the  catalyst  molecule 
plane,  the  other  O  atom  is  farther  away  from  the  metal  atom, 
and  the  projection  of  the  0—0  bond  axis  on  the  catalyst 
molecule  plane  is  aligned  with  the  bisecting  line  of  the  two 
neighboring  bonds  between  the  central  metal  atom  and  the  two 
nearest  N  atoms.  In  the  most  stable  adsorption  of  H2O,  the  O 
atom  lies  right  above  the  central  metal  atom  over  the  catalyst 
molecule  plane,  the  H20  molecule  is  nearly  parallel  to  the 
catalyst  molecule  plane,  and  the  projection  of  the  H— O— H  bond 
angle  on  the  catalyst  molecule  plane  is  bisected  by  the 
connecting  line  between  the  central  metal  atom  and  one  of  its 
nearest  N  atoms.  In  the  most  stable  adsorption  of  OH,  the  O 
atom  lies  right  above  the  central  metal  atom  over  the  catalyst 
molecule  plane,  the  H  atom  is  extended  away  from  the  metal 
atom,  and  the  projection  of  the  O— H  bond  axis  on  the  catalyst 
molecule  plane  is  aligned  with  the  bond  between  the  central 
metal  atom  and  one  of  its  nearest  N  atoms. 

Furthermore,  we  optimized  all  the  possible  structures  of 
HOOH  and  H200  molecules  adsorbed  on  the  catalyst  molecules 
(shown  in  Figure  8).  We  determined  the  adsorption  energy  for 
HOOH  on  FePc  to  be  —0.99  (Figure  8a),  —1.01  (Figure  8b), 
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Figure  6.  The  performance  recovery  test  for  FePc  preadsorbed  on  a 
GC  rotating  disk  electrode  surface  at  a  rotating  rate  of  2500  rpm  in 
oxygen-saturated  0. 1  M  NaOH  at  room  temperature  with  a  scanning 
rate  10  mV  s_1.  In  the  figure,  the  potentials  are  given  with  reference  to 
both  Hg/HgO  electrode  and  normal  hydrogen  electrode  (NHE).  Note 
that  the  computed  potential  of  the  Hg/HgO  electrode  versus  an  NHE 
at  25  °C  and  1  atm  is  0.164  V  for  0.1  M  NaOH  solution. 

TABLE  2:  Calculated  Adsorption  Energy  (in  units  of  eV)  of 
Reactant  (O2  and  H20),  Product  (OH),  and  Reaction 
Intermediate  (HOOH  and  H200)  Molecules  Adsorbed  on 
FePc  and  CoPc  Catalyst  Molecules  Using  the  First-Principle 
DFT  Method" 

02  H20  OH  HOOH  H200 

FePc  -1.16  -1.05  -3.41  -1.05  -1.03 

CoPc  -0.40  -0.32  -2.36  -0.24  -0.47 

Pt  (111)  -0.79  -0.20  -2.39 

a  For  comparison,  the  adsorption  energy  of  02,  H20,  and  OH 
molecules  on  a  Pt  (111)  surface  is  also  listed.  These  adsorption 
energies  on  the  Pt  (111)  surface  were  calculated  with  VASP 
software  and  the  plane  wave  pseudopotential  method. 

—  1.05  (Figure  8c),  and  —0.32  eV  (Figure  8d).  We  found  the 
adsorption  energy  for  H200  on  FePc  to  be  —0.94  (Figure  8e) 
and  —1.03  eV  (Figure  8f).  Thus,  the  configuration  shown  in 
Figure  8c  is  the  lowest  energy  structure  for  H202  adsorbed  on 
the  FePc  molecule.  For  HOOH  on  CoPc,  we  found  that  parts  a 
and  b  of  Figure  8  are  energy  favorable  configurations,  but  parts 
c  and  d  of  Figure  8  are  not.  The  adsorption  energy  for  HOOH 
on  CoPc  was  predicted  to  be  —0.24  (Figure  8a)  and  —0.22  eV 
(Figure  8b).  We  calculated  the  adsorption  energies  of  H200 
on  CoPc  to  be  —0.47  (Figure  8e)  and  —0.45  eV  (Figure  8f). 
Consequently,  the  configuration  shown  in  Figure  8e  turns  out 
to  be  the  lowest  energy  structure  for  H202  adsorbed  on  the  CoPc 
molecule. 

IV.  Discussion 

Our  theoretical  calculation  results  reveal  molecular  details 
about  the  ORR  mechanisms  on  FePc  and  CoPc  catalyst 
molecules,  which  enables  us  to  better  interpret  our  experimental 
observations,  such  as  the  ORR  onset  potentials,  reaction 
pathways  (4e~  or  2e~),  and  stability  of  these  catalysts. 

IV.l.  Oxygen  Adsorption  To  Be  the  Key  Step  To  Control 
the  ORR  Onset  Potential.  First,  we  compare  the  adsorption 
energy  of  the  two  reactant  (02  and  H20)  molecules  on  the 
surface  of  the  three  catalysts  (FePc,  CoPc,  and  Pt  (111)).  Our 
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(C)  (d) 

Figure  7.  Optimized  lowest  energy  structure  of  (a)  isolated  FePc  or 
CoPc  catalyst  molecules,  (b)  Ot  adsorbed  on  FePc  or  CoPc,  (c)  H2O 
adsorbed  on  FePc  or  CoPc,  and  (d)  OF!  adsorbed  on  FePc  or  CoPc. 
The  central  yellow  ball  represents  a  metal  Fe  or  Co  atom,  blue  balls 
represent  N  atoms,  gray  balls  represent  C  atoms,  and  light  white  balls 
represent  H  atoms. 


results  in  Table  2  point  out  that  the  O2  molecule  would  have  a 
lower  adsorption  energy,  and  hence,  it  would  bind  more  strongly 
to  the  three  catalysts  than  the  H20  molecule.  Since  FePc  and 
CoPc  molecules  only  provide  one  adsorption  site  (the  central 
metal  atom),  this  result  suggests  that  the  adsorption  of  the  O2 
(not  the  H20)  molecule  on  FePc  and  CoPc  would  be  the  first 
step  of  the  ORR  in  an  alkaline  solution.  Moreover,  our  results 
in  Table  2  show  that  the  adsorption  energies  of  O2  on  the 
catalysts  follow  the  increasing  order  of  FePc,  Pt  ( 1 1 1),  and  CoPc, 
which  correlates  well  with  our  experimental  observations  (shown 
in  Figures  1  and  3).  FePc/C  catalyst  shows  the  most  positive 
ORR  onset  potential,  which  is  more  than  the  Pt/C  catalyst  and 
much  more  than  the  CoPc/C  catalyst.  On  the  basis  of  the  slopes 
of  the  Tafel  plots  (Figure  4),  the  rate-determining  step  for  ORR 
was  considered  to  be  similar  on  both  the  FePc/C  and  the  Pt/C 
electrodes.  It  is  well  believed  that  O2  adsorption  is  a  limiting 
step  for  the  ORR  on  Pt  surface.  Therefore,  we  conclude  that 
the  O2  adsorption  process  is  a  limiting  step  affecting  the  kinetics 
of  the  ORR  on  FePc  and  CoPc:  the  lower  the  O2  adsorption 
energy,  the  higher  the  kinetics  of  the  ORR  on  the  molecular 
catalysts.  Yeager  et  al.37  had  postulated  earlier  that  the  affinity 
of  Fe(II)Pc  for  O2  may  be  very  high  compared  to  that  for  the 
reaction  intermediates,  products,  or  other  solution-phase  com¬ 
ponents  and  suggested  that  this  strong  interaction  with  O2 
compared  to  that  with  other  species  can  be  of  critical  importance 
to  the  electrocatalytic  activity  involving  a  single  type  of  iron 
site.  Consequently,  our  finding  in  this  work  supports  such  a 
proposition  that  O2  adsorption  plays  a  key  role  in  the  activation 
and  reduction  of  oxygen. 

IV.2.  H2O2  Adsorption  To  Distinguish  the  4e  or  2e 
Routes  of  the  ORR  on  FePc  vs.  CoPc.  The  ORR  would  assume 
the  4e~  route  on  FePc  molecule  catalysts  but  take  the  2e~  routes 
on  CoPc  molecule  catalysts.  As  shown  in  Figure  3,  we 


(e)  (f) 

Figure  8.  Optimized  structural  configurations  for  HOOH  (a— d)  and 
FFOO  (e  and  f)  molecules  adsorbed  on  FePc/CoPc  molecules.  The 
central  yellow  ball  represents  the  metal  Fe/Co  atom,  the  two  red  balls 
represent  O  atoms,  blue  balls  represent  N  atoms,  gray  balls  represent 
C  atoms,  and  light  white  balls  represent  H  atoms. 

demonstrated  that  no  H2O2  was  detected  on  the  FePc/C 
electrodes,  but  significant  amounts  of  FFCF  were  observed  on 
the  CoPc/C  electrodes  with  the  starting  potential  the  same  as 
the  ORR  onset  potential  of  the  disk  electrode.  Our  theoretical 
calculations  provide  an  explanation  as  to  why  this  happens. 

On  FePc  molecules  at  the  lowest  energy  state,  H202  will  have 
a  geometric  structure  as  shown  in  Figure  8c.  DFT  calculation 
predicts  in  Figure  8c  that  the  two  O  atoms  in  an  H2O2  molecule 
bind  to  the  common  central  Fe  atom  with  a  distance  of  1.820 
A,  bind  to  individual  FI  atoms  with  a  distance  of  0.972  A,  but 
are  separated  from  each  other  by  a  distance  of  2.368  A.  In  an 
isolated  HOOH  molecule,  the  0—0  bond  length  is  calculated 
to  be  1.487  A.  Hence,  0—0  bond  breaking,  which  is  a 
prerequisite  for  four-electron  ORRs,  will  occur  during  the 
process  of  H2O2  adsorption  on  FePc  molecules.  In  contrast. 
Figure  8c  shows  an  unfavorite  adsorption  configuration  for 
HOOH  on  CoPc  molecules.  Instead,  the  H202  molecule  prefers 
to  have  the  configuration  of  Figure  8e  when  adsorbed  on  CoPc. 
In  that  optimized  configuration,  the  two  O  atoms  in  the  FFO? 
molecule  are  only  distanced  by  1.536  A,  which  is  very  close  to 
the  equilibrium  0—0  bond  length  of  1.501  A  in  an  isolated 
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H200  molecule.  Clearly,  0—0  bond  breaking  does  not  happen 
during  the  process  of  H202  adsorption  on  CoPc  molecules. 
Therefore,  our  DFT  calculations  about  the  H202  adsorption 
process  strongly  indicate  that  the  breaking  of  the  0—0  bond 
accounts  for  the  4e  ORR  on  the  FePc  catalyst  molecule,  while 
the  lack  of  the  0—0  bond  breaking  leads  to  the  2c  ORR  on 
the  CoPc  catalyst  molecule. 

IV.3.  OH  Adsorption  To  Be  a  Key  Factor  for  Cycling 
Stability  of  FePc  and  CoPc.  Compared  to  02,  OH  will  have 
even  lower  adsorption  energy  and  hence  a  more  energetically 
favorable  adsorption  on  FePc  and  CoPc  molecules  based  on 
our  DFT  calculation  results  given  in  Table  2.  Since  radical  OH 
is  a  product  of  the  ORR  regardless  of  the  reaction  paths  (4c 
or  2e  ),  the  generated  OH  might  occupy  those  active  reaction 
sites  on  the  surface  of  the  catalysts  after  the  ORR  and  prevent 
further  adsorption  of  02.  Conceivably,  this  will  cause  a 
degradation  of  the  performance  of  the  molecule  catalysts.  Using 
the  data  in  Table  2,  we  derived  that  the  difference  in  the 
adsorption  energy  of  02  and  OH  is  about  2.25  eV  on  FePc  and 
1.96  eV  on  CoPc.  Thus,  this  implies  that  the  performance 
degradation  caused  by  OH  would  be  more  severe  for  FePc  than 
it  would  be  for  CoPc  molecule  catalysts.  In  Figure  5,  we  showed 
that  ORR  activity  of  FePc  declined  much  more  rapidly  than 
that  on  CoPc  with  either  continually  cycling  or  holding  at  a 
constant  potential.  Thus,  experimental  observation  on  the 
stability  of  catalysts  is  consistent  with  what  was  predicted  by 
the  DFT  calculations  based  on  OH  adsorption.  Furthermore,  we 
showed  in  Figure  6  that  the  catalytic  activity  of  the  FePc/C  can 
be  partially  recovered  simply  by  leaving  the  catalysts  in  a  0. 1 
M  NaOH  solution  for  some  time.  It  is  reasonable  to  interpret 
the  partial  recovery  of  the  catalytic  activity  as  a  result  of  the 
removal  of  the  adsorbed  OH  on  FePc  catalysts  owing  to  OH 
diffusion  processes  in  alkaline  solution.  After  the  adsorbed  OH 
is  diffused  away  from  the  FePc  molecules,  the  catalysts  regain 
their  active  sites  for  oxygen  molecule  adsorption.  Thus,  the 
results  in  Figure  6  would  further  support  our  proposal  that  the 
strong  OH  adsorption  on  the  FePc  molecule  is  responsible  for 
the  performance  degradation  of  the  FePc  catalyst.  Moreover, 
the  fact  that  the  catalytic  performance  could  not  be  fully 
recovered  implies  that  some  irreversible  processes,  like  molecule 
decomposing  and/or  central  transition  metal  atom  removal, 
would  also  occur  during  the  cycle  tests. 

V.  Conclusions 

Using  RDE  and  RRDE  techniques  and  DFT  calculations,  we 
elucidated  the  mechanisms  of  the  ORR  on  carbon-supported 
FePc  and  CoPc  catalysts.  In  our  electrochemical  measurements, 
FePc/C  displays  not  only  more  than  100  mV  positive  ORR  onset 
potential  and  half-wave  potential,  but  also  more  than  double 
the  limiting  current  than  CoPc/C.  Our  DFT  calculations  gave 
02  adsorption  energy  of  —1.16  eV  on  FePc  and  —0.40  eV  on 
CoPc  molecules.  Correlating  the  experiment  with  theory,  we 
infer  that  oxygen  adsorption  on  the  catalysts  is  the  key  step  to 
control  the  ORR  onset  potential.  For  the  molecular  catalysts, 
the  lower  the  02  adsorption  energy,  the  higher  the  kinetics  of 
the  ORR  that  can  be  expected.  Moreover,  our  electrochemical 
measurements  show  that  FePc/C  catalysts  were  found  to  promote 
the  desired  4c  ORR  in  alkaline  solutions  while  CoPc/C  only 
facilitates  2c  ORR.  Corresponding  to  that  observation,  our  DFT 
calculations  indicate  that  the  ORR  pathways,  2c  or  4c  ,  are 
mostly  determined  by  the  structures  of  H202  adsorption  on  FePc 
and  CoPc  catalysts.  The  breaking  of  the  0—0  bond  during  H202 
adsorption  accounts  for  the  4e  ORR  on  FePc  catalyst  molecule, 
while  the  lack  of  the  0—0  bond  breaking  process  leads  to  the 


2c  ORR  on  CoPc  catalyst  molecule.  Third,  we  believe  that 
OH  adsorption  on  the  catalysts  is  an  important  factor  to 
determine  the  stability  of  the  FePc  and  CoPc  catalysts.  It  was 
predicted  from  DFT  calculations  that  OH  adsorption  is  more 
energetically  favorable  than  02  adsorption  on  both  FePc  and 
CoPc  catalyst  molecules.  This  will  cause  the  produced  OH 
molecule  to  occupy  the  active  sites  of  catalysts  and  hamper 
further  progress  (02  adsorption)  of  ORR.  Our  DFT  result  implies 
that  this  degradation  process  should  be  more  severe  for  FePc 
than  CoPc,  due  to  a  larger  energy  gain  when  adsorbing  OH 
instead  of  02  on  the  FePc  surface.  Concurring  with  DFT  work, 
our  experimental  stability  tests  showed  clearly  that  the  FePc/C 
catalyst  will  lose  its  performance  much  more  quickly  than  the 
CoPc/C  catalyst. 

Thus,  we  have  demonstrated  that  the  combination  of  elec¬ 
trochemical  measurements  (RDE  and  RRDE)  and  theoretical 
modeling  is  a  powerful  approach  in  identifying  high-perfor¬ 
mance  catalysts  for  ORR  and  investigating  the  ORR  mecha¬ 
nisms.  Our  study  suggests  that  the  FePc/C  is  a  promising 
nonprecious  electrocatalyst  and  can  be  applied  in  alkaline  fuel 
cells.  The  ORR  on  FePc/C  catalyst  has  a  quite  positive  onset 
potential  and  half-wave  potential,  and  adopts  a  high  energy 
capacity  4c  reaction  path.  However,  the  FePc/C  catalyst  will 
lose  its  catalytic  activity  rapidly  under  fuel  cell  working 
conditions.  Understanding  the  mechanisms  underlying  those 
observations  as  presented  in  this  paper  is  essential  for  us  to 
further  tailor  the  chemistry  of  the  catalyst  molecules  for  the 
purpose  of  finding  an  economic,  active,  and  durable  catalyst 
for  SAFCs. 
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